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Abstract 

Abelian fermion models described by the SLAC action are considered on a finite 
2d lattice. It is shown that modification of these models by introducing additional 
Pauli - Villars regularization supresses nonlocal effects and provides agreement 
with the continuum results in vectorial U(l) models. In the case of chiral fermions 
the phase of the determinant differs from the continuum one. 
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1 Introduction 



In a recent paper [[Xj] we considered lattice fermion models on the 2d lattice 
with Wilson action improved by Pauli - Villars (PV) regularization. It was 
shown that in spite of the chiral symmetry breaking for finite lattice spacing 
a, its effects are supressed by PV regularization and the model provides a 
good agreement with the continuum results both in the perturbative and 
nonperturbative region. Nevertheless lack of chiral symmetry for finite a is 
not quite harmless as to get a good agreement for large external fields one 
needs big lattices. 

So it would be highly desirable to have a formulation which preserves 
chiral symmetry for a finite lattice spacing as well. 

"No-go" theorem j2|] forbids any local formulation preserving chiral sym- 
metry, therefore we consider a nonlocal SLAC model proposed originally in 
paper [j3j]. It is known that the original formulation of the model requires 
introduction of nonlocal counterterms [[3], |5j, which makes it unacceptable 
for practical calculations. However the model can be improved in the same 
way as it has been done for Wilson action by introducing additional gauge 
invariant PV regularization which supresses the contribution of momenta 
close to the edge of the Brillouin zone. It was shown || that for anomaly 
free models on an infinite lattice in the framework of perturbation theory 
all nonlocal effects can be supressed in this way and one gets a manifestly 
chiral invariant formulation of anomaly free models. 

The propose of this paper is to check these results both perturbatively 
and nonperturbatively for the models on the 2d finite lattice. We found that 
in the case of vectorial models this approach works leading to a reasonable 
agreement with continuum results. However for chiral fermions, even when 
perturbative anomaly is absent, we observed a discrepancy in the value 
of the phase of continuum and lattice determinants. The origin of this 
discrepancy is discussed. 
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2 Vectorial lattice model 



In this section we consider the vectorial model described by the action 

_ y» 
Ivs= E^W7#"^xp[i E A »( z M(y)- (!) 

x,y,/x z^ — 

Here —N/2 + 1 < x^ < N/2, n = 0, 1. The lattice spacing is chosen to be 
equal to 1. We restrict ourselves by the case of a constant external fields 

= TV M = const " 

The Fermi field ^ satisfies antiperiodic boundary conditions. V^x) is the 
lattice (SLAC) derivative: 



1 ^ ww N 27Ti/ 1 



= t^2 E_ i^m(p) ex P -ft- [P 

where 



= lOrfa ~ l\ ~ N / 2 + 1 < *V " ™N < N/2, m = 0, ±1, 



N v f 2 

The action (g) is gauge invariant but it is not local. 

A straightforward calculation gives for the vectorial fermion determinant 
normalized to 1 at h = 0, the following expression: 

where B 2 (p, h) = £ -B^(p, /i), B M is the Fourier transform of the covariant 
derivative 

*,(p.fc) = £ e (-ir +l Sinf (3) 

iV *=-JV/2+l sm j\T 

The dependence of B^/2tt on (p M — — l/2)/N is presented at Fig.l. One 
sees that on a finite lattice, the saw-tooth form of the SLAC derivative is 
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Fig.l. Covariant derivative B(p, K)/2tt as function of (p — h — 1/2) /N 
at TV = 16 and 160 



smoothen and contrary to the infinite lattice case may be approximated by 
a continuous curve. This fact will be important for further discussion. 

The lattice determinant (|2J) has to be compared with the known result 
for the continuum theory [j7], @, H 

D vc = e- 2 ^ n \F[n,h]F[n, -h}\\ (4) 

n=l 

where 

^ e -27r(n-l/2)+27ri(/i +ifti) 
F t n ' h \ = I ~ e -27r(n-l/2) • 

The determinants (@), (f|) satisfy the following symmetry properties: 

D[h , hi] = D[h u h ] = D[-h , hi] = D[h + n , hi + m], 

where no, ni = 0, ±1, ±2, . . .. Therefore it is sufficient to consider the fields 
h^ only in the range < ho < hi < 1/2. 
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Fig. 2. Vectorial diagrams as functions of N at external momentum k = 0: 
1-U VS {0); 2, 3 - with four lines: 2-r 000 o(0), 3 - r ioi(0) 



Using the eq. (@) one can easily write the expression for the diagrams 
with two and more external /i-lines. It looks as follows: 



2 N £ I Bg(p,h) + B < ,(p,ft)B>,h) 

dhydhy p =-n/2+i 



W^p, h)B'^ h)B„(p, h)B' v { P , h) x 

rn,/ r, } 2 ( 5 ) 



[B*(p,h)]' 

where ' means a derivative with respect to h. We wish to study the be- 
haviour of different diagrams at N — ► oo. Computer simulations of the 
second and fourth order diagrams as functions of TV were performed using 
eqs (|5|) and ([3D. The results are presented at Fig. 2. One sees that the polar- 
ization operator U V s(0) and the diagram with four external lines r ooo(0) 
diverge like N. These calculations show that the SLAC action does not 
provide the correct continuum limit. A mass renormalization which one 



5 



could expect on the basis of power counting in the continuum theory does 
not save the situation as the diagrams with more than two external lines 
diverge in the limit N— >oo. 

We shall try to improve the model by introducing additional PV regu- 
larization according to ||. The regularized action looks as follows: 



Here 4> r are Bose and Fermi PV fields having the same spinorial and internal 
structure as ip. 

The estimates of asymptotic behaviour of different diagrams, analogous 
to the ones presented above, show that if one uses one PV field the dia- 
grams have asymptotics ~M 2 7V" 2 . Therefore to supress the contribution of 
momenta close to the edge of the Brillouin zone which are responsible for 
nonlocal effects, one has to choose M <C jj. Such small values of PV field 
masses are not acceptable as they are comparable with masses of physical 
particles. Moreover, as we shall see from numerical analysis, the model with 
one PV field is very sensitive to the particular choice of M and therefore 
the results are not stable. 

To get reliable results one needs to introduce at least three PV fields. 
In this case the diagrams have asymptotic behaviour M A N 2 , and choosing 

<C M <C -7= one can supress the contribution of momenta close to the 
edge of the Brillouin zone. 

Below we present the results of numerical calculations for the cases of 1 
and 3 PV fields. 

The regularized determinant has a form 



Ivr — Ivs + Ipv 



ipv = E E ^rOh/*£VO - 2/) ex p[i E A n( z )](f>r(y)+ 




x 



D VR = D V s{h}D PV [h}, 
where Dvs[h] is given by eq. (j2j) and Dpy is defined as follows: 



(6) 
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Here c r = 1 or —1 corresponds to the case of Fermi or Bose PV field. 
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Fig. 3. Vectorial determinants Dy as functions of M at ho = h\ = 0.2: 

1 - Dye on the torus; 2, 3, 4 - Dvr with PV fields: 

2 - 7V = 32, 3 PV fields; 3, 4 - 7V = 160: 

3 - 1 PV field, 4 - 3 PV fields 

The results of calculations are presented at Fig. 3. One sees that in the 
case of one PV fields the agreement is achieved for iV=160 only at one 
particular value of Mo=0.026 - 0.033. Even small variation in the value of 
M leads to a large discrepancy between lattice and continuum results. In 
the case of 3 PV fields the lattice results agree with the continuum in the 
rather big interval of values of PV fields masses. For 7V=160 the agreement 
is observed for Mi = 0.04 < M < = 0.2. These values are practically 
independent of h, except for the case — > 0.5. Comparing the data at 
iV=32 and 7V=160 one notes that if N grows the mass values decrease and 
the minimal value of M\ at 3 PV fields behaves like iV -3 / 4 . 

Let us discuss the behaviour of the lattice determinant in the case of 
3 PV fields in some more detailes. Due to the fact that the covariant 
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derivative B^(p, h) sharply goes to zero at the border of the Brillouin zone, 
it is instructive to separate the products in eqs (@), (0) into three parts: 

D VR = D in D b D a . (7) 

Here D m is a product over p excluding external points, is the product 
over p belonging to the edges of the Brillouin zone and D a is the product 
of p corresponding to the vertices of the zone. 

Let us start with D m . According to the behaviour of B^(p, h) (see Fig.l), 
in this region difference between regularized diagrams and the continuous 
ones is of order 1/MN. Near the border of this region B^(p,h) and its 
derivatives are of order 1. It allows to expand the corresponding terms in 
regularized diagrams over M 2 . In the case of 3 PV fields the first nonvan- 
ishing term is ~M 4 leading to the asimptotic behaviour ~M 4 N 2 . Summa- 
rizing these estimates we get 

D VR = D vc (1 + (1/MN) + 0(M 4 7V 2 )) , N -> oo, (8) 

when £ <C M < 4-. 



For the further estimates of values D a and D^, using numerical data 
presented at Fig.l, we approximate the covariant derivative B^p, h) by the 
following two lines (ignoring oscillations): 

2ir 

B,(P, h) « —{p, - h, - 1/2), PfX = -N/2 + 2, . . . , N/2 - 1; 

B^p + N/2, h) « -2n(p, - hp - 1/2), p fl = 0, 1. (9) 

The value of D a at | \ < 1/2 has the form following from the eqs (§) and 
®: 

D a = D a0 [h , hi]D a0 [h , -h^Daolhu h ]D a0 [h h -h ]/D 4 a0 [0, 0]. 

Here 

JJaO[tlO, fll\ — ~ 



[(1 p) 2 +(M/27T)2 
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2 2 

where (| — p) = E (| — J?^) . One sees that at all values of h in the interval 
\hn\ < 1/2 except the trivial case |/io| = = 1/2 when Dvr = Dvc = ® 



D a ^l atM(JV)->0. 



(10) 



Now we consider the value D b . Its expression at \h^\ < 1/2 looks as 
follows: 

D b = D b0 [h 0j h x ] D b0 [h , -h^Dwlhu h }D b0 [h u -/i ]/£> 6 4 [0, 0]. 



Here 



L J Po =-^/2 + 2 [G[p ,/ i ] + (M/27T)2] 2 



where G[po, ^] = (po ~ ho - §) + (§ - foi] 

Let us extend in the last equation from the interval —N/2 + 2 < £> < 
TV/ 2 — 1 to — oo <po < oo. In additional domain ]?o — \ > 
behaves like 



iV-1 



the value D 



1 Kiv-i 

PO-2 



--N PO 



Using this estimation one can transform D b Q to the form 



where M] = ch27riVy (± - /h) + (ff) - cos2tt (/i , 2 



2 . / m \ 2 / , | 

2, 

According to this formula, we get for the value D b at | — ^> M the 
following expression: 

M 4 7V 



= exp 



32tt 3 



(1/2 - /n)~ 3 + (1/2 + /h)- 3 + 



(1/2 - h y 3 + (1/2 + 



,-3 



32 



(11) 



It follows from eqs (0), (@), ([TDD, (03) that for ^< M<c 7/v our regular- 
ized lattice model agrees with the continuum toron model for the fields 
in the interval \\ ± h^\ >(M 4 7V) 1/3 . 
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Fig. 4. Vectorial determinants Dy as functions of h\ at ho = 0.2: 

1 - Dye on the torus; 2, 3, 4, 5 - £Vi? with 3 PV fields: 
2, 3-7V = 32: 2-M = Mi = 0.15, 3-M = M 2 = 0.4; 
4, 5-7V=160: 4-M = Mi = 0.04, 5-M = M 2 = 0.15 



In the region (M^N) 1 ^ 3 ^> ||±/i M | ^> M, the lattice determinant decreases 
sharply when \h^\ — > 1/2 and the agreement with the continuum model is 
lost. This effect is due to vanishing of the covariant derivative B fi (p, h) at 
the border of the Brillouin zone, and as our calculations show it comes from 
the contribution of the edges of the Brillouin zone. 

If we denote the value of h at which Dy^ and Dye start to differ by more 
than 5% by h*, i.e when DvrJ Dyc = £,=0.9b, it follows from eq. that 

1 /il 2 V327r 3 ln^-V ' 
One sees that in the limit N — > oo the regularized model agrees with the 
continuum toron model in the whole interval of h values. 

These analytic results are in a good agreement with numerical calcula- 
tions of DyR and Dye according to eqs (|6j), ([§). At Fig. 4 the values of Dyn 
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and Dye calculated for different h and M are presented. One sees that 
for 7V=160 3 PV fields provide a good agreement with the continuum for 
\h\<0A. 

3 Phase of the lattice chiral determinant 

In this section we consider the possibility to use the SLAC action for 
regularizing chiral U(l) model. Obviously such a regularization could be 
successful only for anomaly free models, as SLAC action is chiral invariant 
and therefore cannot reproduce anomaly. However our discussion in the 
previous section shows that even in the case when perturbative anomalies 
are compensated as in 11112 model, this regularization most probably fails. 
Indeed the possibility to supress the contribution of momenta close to the 
border of the Brillouin zone in anomaly free models by introducing gauge 
invariant vectorial interaction of PV fields is related to the fact that usually 
we deal with a finite number of divergent diagrams, whose sum is anomaly 
free models is parity conserving. For example in the 2d U(l) model with 4 
positive chirality fermions with charge 1 and one negative chirality fermion 
with charge 2 the sum of anomalous second order diagrams is purely vector 
like. For nonzero external momenta all higher order diagrams are conver- 
gent and the contribution to them of momenta close to n/a is negligible. 
This is the reason why PV regularization was successfully applied to the 
SLAC action on the infinite lattice, see [j6|]. However in the toron model with 
SLAC action, as was shown in the previous section, there are "divergent" 
diagrams with more than two external lines. That means the contribution of 
boundary momenta ~ - is important not only for two-point diagram, but 
for other diagrams as well. The sum of the diagrams with more than 2 ex- 
ternal lines is not vector-like and therefore cannot be supressed by vectorial 
PV interaction. Using PV regularization one can easily obtain the agree- 
ment with the continuum case for the modulus of determinant. However 
vectorial PV interaction does not influence the phase of the determinant. 
Therefore the SLAC action can provide a correct phase only if for some 
reasons the contribution of boundary momenta cancel by ifself. Moreover it 
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is easy to show that nonzero contribution to phase comes from the diagrams 
with more than two external lines and hence anomaly cancelation does not 
help. 

So it is sufficient to check the phase of a model with one positive chiral- 
ity fermion described by the standard SLAC action ([I]), where one has to 
consider ip as the Weyl fermion 

No PV fields are needed. 

A straightforward calculation gives the result 

+S p ^ 1 /2+1 Bob,0)+i5 1 (p,0)' 1 } 

where B^(p, h) is the covariant derivative defined by the eq. (|3|) (see also 
Fig.l). 

The corresponding expression for the determinant in continuum theory 
on the torus looks as follows J7|, @, |9j]: 

oo 

D +c = e [7rh ^ ha+ih ^ J] F[n, h]F[n, -h], (13) 

n=l 

where F[n,h] is defined by eq (f|). 

The determinants (|T2|), Q13Q satisfy the following symmetry properties: 

D + [h , /ii] = D* + [h u ho] = D* + [-h , /ij. 

Due to chiral invariance of the SLAC action the lattice determinant ( |I3f) is 
periodic in /i: 

D +s [ho, hi] = D +S [h + n , /ii + nj, n , n\ = 0, ±1, ±2, . . . . 

The continuum theory is anomalous and the corresponding determinant 
satisfies the condition 

D +C [h + n , h + ni ] = ^ n ^- n ^D +c [h , h]. 
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It follows that one can hope at most on the agreement only for the fields 
hfj, in the interval |/i M |<0.5. It is sufficient to consider < ho < hi <0.5. 

Let us firstly make analytic estimates of Arg D + $. In the previous section 
we showed that due to sharp decreasing of B fi (p, h) near the border of 
the Brillouin zone this region gives a considerable contribution which may 
spoil the agreement with the continuum case. For that reason we shall 
study the behaviour of Arg D + $ i n the interier and boundary part of the 
zone separately. Using the notations introduced in the previous section we 
present Arg D + s as a sum 

Arg D + s = Arg D +in + Arg D +b + Arg D +a mod 2tt. (14) 

To estimate the separate terms we adopt a linear approximation (0) of the 
covariant derivative B^(p, h). Eq. ( |12|) leads to the following expression for 
different terms: 

ArgD +in = Yl arctg — f f ^ $4- mod 2tt, (15) 

p=-N/2+2 \ [Pfi - n v - 2) \Pf* ~ 2) 



Arg D +a = arctg 



4 {hi - h\) - lQhlh\ + 1 

AigD +b = Arg D +b0 [h , hi] + ArgD +b0 [-h , -hi] + 
+ AigD +io [-hi,ho] + Arg D +b0 [h 1: -h ] mod 27r, 



n N/ ^ 2hoN(hi-\, 

Arg D +b0 [h , hi] = E arctg- -2 — — mod 2tt. 

po=i (po -5) + A^ 2 (/ii - i) - ^ 

Due to the symmetry properties of the determinant D + $ only the dia- 
grams with more than 2 external lines contribute to Arg D +m . These lattice 
diagrams differ from the corresponding continuum diagrams by the terms 
of order 1/N. Therefore for | /i^ | < 0.5 we have 

ArgL> +in ^ ArgL> +c , N^oo. (16) 

To estimate Arg D +b o[ho, hi] we expand the arctg a; in the Taylor series. 
Keeping only the first term we replace the sum by the integral and substitute 
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Fig. 5. Arguments of positive chiral determinants ArgD + 
as functions of hi at ho = QA: 

1 - Arg D + c on the torus; 2, 3 - Arg D + s on the lattice: 

2 - computed by eq. ( |T2|) at A^ = 32 and 160, 

3 - estimated by eqs ([[4]) ~ (O) at Af ^ oo 



this expresion into the formula for ArgD+^o. In this way we get 



Arg D + b = 2hi arctg 



2h f 



l-2hl 



2ho arctg 



2h x 



1-2/if 

The calculations show that for all h in the interval under consideration 



(17) 



\ArgD +b \ <C | Arg D +in \, N^oo. 
Finally, calculating Arg D +a with the help of eq. ([TBI) we find 

ArgL> +a = ArgL> +in , N^oo. 



Our estimates show that the contribution of the edges of the Brillouin 
zone to the Arg D + $ is small, but the contribution of the vertices is approx- 
imately equal to the contribution of the interier and their sum is two times 
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bigger than the continuum value: 



AigD + s = 2kigD +c . 



So we have here some kind of doubling phenomenon which is somewhat 
reminiscent to the phenomenon observed by Bodwin and Kovacs [jiTlfl in 
Rabin's formulation of Schwinger model 



Analytic estimates given above are in a good agreement with the results 
of numerical calculations presented at Fig. 5. For all values of h Arg D + s is 
two times bigger than the corresponding continuum value. 



4 Discussion 

In this paper we analyzed the U(l) model on a finite 2d lattice described 
by the gauge invariant nonlocal SLAC action. It was shown that in the case 
of a constant gauge field (toron model) SLAC action generates infinite series 
of divergent (in the limit N — > oo) diagrams with more than two external 
lines and the value of the lattice fermion determinant does not agree with 
the known exact result for the continuum toron model. Modification of 
the SLAC action by introducing additional PV regularization cures this 
decease in the case of vectorial interaction. A minimal number of PV fields 
necessary for such regularization is equal to three. 

The peculiar feature of the SLAC action on a finite lattice is large contri- 
bution of momenta corresponding to the vertices of the Brillouin zone. In 
the case of the vectorial interaction contribution of these momenta results 
in sharp decreasing of the lattice determinant in the narrow region near 
\hfj, \ — 1/2. However in the limit N — > oo the width of this region tends to 
zero and the model agrees with the continuum result for all \h^\ < 1/2. 

In the case of axial interaction the existence of divergent diagrams with 
more than two external lines does not allow to cancel the contribution of 
momenta \p\ ~^ by introducing PV regularization. As in the vectorial case 
the boundary momenta give a large anomalous contribution. As a result 
the phase of the lattice chiral determinant is two times bigger than the 
corresponding continuum phase. 



15 



Our results show that SLAC action supplemented by additional PV reg- 
ularization solves successfully the problem of fermion spectrum doubling in 
the case of vectorial interaction, but fails to describe correctly chiral models 
interacting with a constant gauge field. 
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